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We have directly measured the energy threshold and efficiency for bubble nucleation from iodine 
recoils in a CF3I bubble chamber in the energy range of interest for a dark matter search. These 
interactions cannot be probed by standard neutron calibration methods, so we develop a new tech- 
nique by observing the elastic scattering of 12 GeV/c negative pions. The pions are tracked with a 
silicon pixel telescope and the reconstructed scattering angle provides a measure of the nuclear recoil 
kinetic energy. The bubble chamber was operated with a nominal threshold of (13.6 ± 0.6) keV. 
Interpretation of the results depends on the response to fluorine and carbon recoils, but in general 
we find agreement with the predictions of the classical bubble nucleation theory. This measurement 
confirms the applicability of CF3I as a target for spin-independent dark matter interactions and 
represents a novel technique for calibration of superheated fluid detectors. 

PACS numbers: 29.40.-n, 95.35. +d, 95.30.Cq, FERMILAB-PUB-10-318-A-CD-E 



Recent years have seen a resurgence in the use of su- 
perheated liquids and bubble chambers as continuously 
sensitive nuclear recoil detectors searching for dark mat- 
ter in the form of Weakly Interacting Massive Parti- 
cles ( WIMPs) [IH5] . At a low degree of superheat, bub- 
ble chambers are insensitive to minimum ionizing back- 
grounds that normally plague WIMP searches but retain 
sensitivity to the nuclear recoils that would be charac- 
teristic of WIMP scattering. In a superheated liquid 
the process of radiation-induced bubble nucleation is de- 
scribed by the classical "hot spike" model |3]. For the 
phase transition to occur, the energy deposited by the 
particle must create a critically sized bubble, requiring 
a minimum energy deposition in a volume smaller than 
the critical bubble. Under mildly superheated conditions, 
the latter requirement renders the bubble chamber insen- 
sitive to minimum ionizing particles. 

The radius of the critical bubble is given by the con- 
dition that the bubble be in (unstable) equilibrium with 
the surrounding superheated fluid [5] . This demands the 
pressure balance 

P b -Pi = ~, (1) 

r c 

where Pb is the pressure inside the bubble, Pi is the pres- 
sure in the liquid, a is the bubble surface tension, and 
r c is the critical bubble radius. The pressure Pj, is fixed 



by the condition that the chemical potential inside and 
outside the bubble be equal, giving 



(A ~ Pi) « (Psat ~ Pi 



Pi - Pv 
Pi 



(2) 



where P sat is the pressure in a saturated system at the 
given temperature, and pi and p v are the liquid and vapor 
densities in the saturated system [BJ. 

The energy threshold for bubble nucleation is the heat 
input necessary to create the critical bubble, given by 



Et = 47rr 2 a 
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Here, T is the temperature of the system, the temper- 
ature derivative is taken at constant chemical potential, 
Pb is the bubble vapor density, and hb and hi are the 
specific enthalpies of the bubble vapor and superheated 
liquid. The three terms give, from left to right, the heat 
necessary to create the bubble surface, the heat needed 
to vaporize the fluid to make the bubble interior, and a 
reversible work term done in expanding the bubble to the 
critical size that must be subtracted to prevent double- 
counting work present in both of the first two terms. To 



good approximation hb — hi may be replaced by the nor- 
mal heat of vaporization at temperature T. The greatest 
uncertainty in determining the thermodynamic Et is the 
relation between the normal surface tension between sat- 
urated liquid and vapor and the surface tension a for 
a very small bubble. This relation is described by Tol- 
man [7j and gives an uncertainty on Et of ~3%. For the 
rest of this paper, we refer to the calculated threshold in 
Eq. (|3| as the Seitz threshold. 

The Seitz model assumes the efficiency for bubble nu- 
cleation is 100% for all interactions that deposit E > Et 
over a volume small compared to the critical bubble. The 
length scales for nuclear recoil cascades in the energy re- 
gion between 5 and 20 keV relevant for a WIMP search 
are similar to the critical radius, so direct calibrations of 
bubble nucleation efficiency are necessary. 

The working fluid discussed in this paper is iodotri- 
fluoromethane or CF3I, which contains two highly sensi- 
tive WIMP target nuclei: fluorine, for spin-dependent 
interactions, and iodine, for spin-independent interac- 
tions. Neutrons are typically used to mimic WIMPs in 
calibrating the nuclear recoil response of a WIMP de- 
tector, and neutron sources have been used to measure 
the nucleation threshold for carbon and fluorine recoils 
in CF3I, CF 3 Br p][8] and C4F10 [9] under various super- 
heat conditions. However, iodine recoils contribute only 
a small fraction to the total neutron-nucleated bubble 
rate in CF3I. Therefore, neutron sources arc an ineffec- 
tive calibration tool for iodine recoils in COUPP. We have 
used heavy daughter nuclei produced in alpha decays as a 
proxy [10], but these are high energy recoils of ^100 keV. 
This paper describes a measurement of bubble nucleation 
efficiency for iodine recoils near our dark matter search 
thresholds. 



EXPERIMENTAL DESIGN 

Our bubble chambers are insensitive to minimum ion- 
izing particles, allowing us to exploit a new calibration 
technique using charged pions as WIMP surrogates to 
produce nuclear recoils by strong elastic scattering. We 
measure the pion scattering angle using silicon pixel de- 
tectors. The nuclear recoil kinetic energy can be calcu- 
lated by E = (p9) 2 /2M on an event by event basis, where 
p is the beam momentum, 9 the scattering angle, and M 
the nuclear mass of the target. 

For a CF3I target, a measured scattering angle corre- 
sponds to a different recoil energy depending on which 
nucleus is involved in the interaction; in this paper, we 
will refer to iodine equivalent recoil energy, Ej e , as the 
energy given to an iodine nucleus for a given pion scat- 
tering angle. For a 12 GeV/c pion beam, approximately 
75% of the rate of pions scattering into angles corre- 
sponding to Ej e between 5 and 20 keV is due to elastic 
scattering on iodine, with smaller contributions from car- 



bon, fluorine, and inelastic scattering [TTj. Therefore, the 
bubble nucleation efficiency for iodine recoils in a bubble 
chamber with Seitz threshold between 5 and 20 keV can 
be inferred from a measurement of the fraction of pion- 
scattering events that nucleate bubbles in the chamber 
as a function of Ei e . 

The measurement was performed in the Fermilab Test 
Beam Facility [T5] using a 12 GeV/c mainly 7r _ beam 
with <Tp/p = 4% and an angular spread of < 1.1 mrad. 
The absolute momentum of the beam is known to 3%. 
The pions were tracked with a silicon pixel telescope [T3] 
consisting of 4 upstream and 4 downstream silicon pixel 
plaquettes, with a spatial coverage of 14 mm x 14 mm. 
The total length of the telescope was 90 cm. The an- 
gular resolution was 0.6 mrad (<r) in the horizontal (x) 
direction and 0.7 mrad in the vertical (y) direction, with 
roughly equal contributions from multiple Coulomb scat- 
tering (MCS) in the target and the spatial resolution of 
the telescope. Plastic scintillators triggered the pixel tele- 
scope on each beam particle. 

A small bubble chamber was designed for this measure- 
ment consisting of a quartz test tube with inner diameter 
10 mm and 1-mm-thick wall, filled with 7 cm 3 of CF 3 I. 
The small size is required to minimize MCS in the short 
radiation length of CF3I (58 mm). The bubble cham- 
ber was operated at a pressure of 30.0 ±0.1 psia and a 
temperature of 34.2 ± 0.2 C with a nominal Seitz thresh- 
old of (13.6 ± 0.6) keV. The iodine equivalent threshold 
scattering angle is 4.7 mrad. An acoustic transducer was 
attached to the top of the test tube to record the acous- 
tic emission produced during bubble formation, provid- 
ing the time of bubble nucleation with ~10 /zs resolu- 
tion. Temperature control was provided by a water bath 
around the bubble chamber. 

Bubble chamber data were taken between March 14 
and March 28, 2012, with a beam flux of ^1000 particles 
per 4-second beam spill with one spill per minute. The 
size of the beam spot was wider than both the bubble 
chamber and the pixel telescope. The chamber was ex- 
panded to the superheated state 22 seconds before the 
arrival of the beam, allowing time for pressure and tem- 
perature transients to dissipate after expansion. The ob- 
servation of bubbles by a 100-Hz video camera system 
created a bubble chamber trigger, causing the video im- 
ages and associated data to be recorded and the chamber 
to be recompressed. After recompression, the chamber 
was dead for the remainder of the beam spill, allowing 
us to collect at most one bubble event per minute. We 
collected about four good single-bubble events per hour, 
with the primary losses due to premature bubble cham- 
ber triggers, bubbles forming outside of the region cov- 
ered by the telescope planes, multiple bubble events and 
large-angle scatters outside the acceptance of the down- 
stream plaquettes. The last two categories are predomi- 
nantly the result of inelastic interactions. Figure [TJ shows 
an example scattering event. 



At the end of the run the CF3I was removed and a 
target empty data set was taken. In addition, data were 
taken in a test run in December 2011 with no target, 
as well as solid targets of quartz, graphite, Teflon or 
(C2F4) n , and crystalline iodine. 



Telescope trigger 



Acoustic signal 




-22.5 -22 -21 .5 

Time (ms) 



o 



o 

Q. 




-400 -200 200 400 

Z position (mm) 



FIG. 1. (Color online) An example event (8 — 6 mrad), in- 
cluding the relative timing of the telescope trigger and acous- 
tic signal, one camera image of the bubble, and the y and 
z positions of the telescope hits. Beam tubes in the water 
bath to minimize the material traversed by the pion beam 
are visible to either side of the bubble chamber 



ANALYSIS 

The primary analysis output is the bubble nucleation 
fraction as a function of Ej e , given by the ratio 



^bub 



N h . 



ub 



(iV tot - AUitO/c 



CF3I 



(4) 



where iVbub is the observed number of pion tracks cre- 
ating single bubbles, N tot is the total number of pion 
tracks, A mu i ti is the number of tracks creating multiple 
bubbles, and /cf3I is the fraction of scatters that occur 
in the active CF3I volume, determined by a comparison 
of the number of scatters in the target-full data set to 
the number in the target-empty data set normalized to 
the number of pion tracks (N emp ): 



/< 



CF3I 
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(5) 



An angular smearing correction is made to N emp to in- 
clude the MCS from the absent CF3I by convolution with 
the standard Gaussian approximation for MCS [T4"] . 

Each pion track is fitted for an upstream and down- 
stream component, with an associated scattering angle 



and 3-D point of closest approach of the two compo- 
nents. The upstream and downstream track segments 
are required to have exactly one hit cluster in at least 
three of the four pixel planes, good fits to straight lines 
(x 2 / v < 4), and to meet in space to within 0.5 mm. To 
exclude pions that passed through little or no CF 3 I, the 
upstream track is required to pass within 4 mm of the 
center of the 10-mm-diameter bubble chamber in the x 
direction. The y location of the track is limited by the 
vertical extent of the pixel planes. Because the uncer- 
tainty on the location of the point of closest approach 
in the beam direction (z) depends strongly on the scat- 
tering angle, we require the z location to be within 3ct z 
of the bubble chamber, where a z is the uncertainty on 
z for each individual event. Events with more than one 
track are rejected. As these track cuts are applied with- 
out regard to the presence of nucleations in the bubble 
chamber, their efficiency applies equally to Ab u b, A to t, 
A^nuiti) and A cmp , and therefore cancels in the final ra- 
tio, r bub . 

The next step is to associate a bubble with a unique 
track using both time and space correlations. The tim- 
ing requirement for correlating a track with a bubble is 
chosen to be 20 < At < 120 /is. The bubble locations 
are reconstructed using standard COUPP techniques [T], 
and the difference between reconstructed bubble position 
and point of closest approach of the track components is 
required to be less than 2.1 mm in the a; direction and less 
than 0.9 mm in the y direction. The combined event ac- 
ceptance of these timing and spatial cuts is 0.958±0.011. 
After these data selection and quality cuts, 350 good sin- 
gle bubble events remain. The final bubble nucleation 
fraction, r DUO , is shown as the points in Fig. [2j 

To disentangle the iodine component from carbon, flu- 
orine and inelastic scattering, we perform a full simu- 
lation using GEANT4.9.5 [IS]- The simulation was vali- 
dated by comparing the simulated scattering angular dis- 
tributions to data for no target, target empty, target full, 
and the solid targets. In all cases, in the MCS-dominated 
small scattering angle region there is good (few per- 
cent) agreement with no adjustable parameters, suggest- 
ing that the telescope geometry is accurately modeled 
in the simulation. In the larger scattering angle region 
dominated by strong elastic scattering, the simulation 
systematically overestimates the observed scattering rate 
by ~ 40%. As this ratio is measured to be the same for 
Teflon (1.45 ± 0.10) and iodine (1.41 ± 0.12), we assume 
that the relative contributions of iodine, fluorine and car- 
bon are being accurately described by the MC. 

A significant systematic uncertainty is introduced by 
our developing understanding of the carbon and fluorine 
recoil nucleation efficiency in this low energy regime. On- 
going studies with ad hoc neutron sources [H] will reduce 
this uncertainty in the future, but here we apply the ex- 
ponential carbon and fluorine efficiency model from |10) : 
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FIG. 2. (Color online) The fraction of pion scattering events 
that produced bubbles as a function of iodine equivalent re- 
coil energy. The solid curves show the simulated contribution 
from individual recoil species (from high to low at 20 keV, 
red for iodine, green for fluorine, and pink for carbon and in- 
elastics), with the blue dashed curve showing the sum. The 
iodine curve shown takes a step function efficiency model for 
iodine recoils using the best fit threshold of (16.8ti'i) keV. 



e(E) = 1 - exp(-a[(E - E T )/E T }), 



(6) 



where E is the nuclear recoil energy and Et is the thresh- 
old energy. In [TO] we found a = 0.15 to be consistent 
with AmBe neutron calibration data. 

We test the hypothesis that the iodine recoil nucle- 
ation efficiency follows the nominal Seitz model of a step 
function (or 3> 1) with 100% efficiency above the Seitz 
threshold by fitting a step function to the data in the 
region 9 < E Ie < 42 keV, allowing E T to float. The fit 
returns Et = (16.8 ^5i) keV, where the error bars are 
statistical. This value is 2.1<r higher than the Seitz model 
threshold (13.6 ± 0.6) keV ± 6%(sys), where the system- 
atic error includes absolute energy scale uncertainties of 
3% in the beam momentum and 1% in the scattering an- 
gle stemming from uncertainty in the z positions of the 
plaquettes. The fit is shown as the dashed blue line in 
Fig. |] 

Figure [3] shows the inferred iodine nucleation efficiency 
as a function of Ej e with the iodine component isolated 
by subtracting the simulated contributions from carbon, 
fluorine and inelastic scatters. The dashed blue curve is 
the best fit step function with E T = (16.8±?;f) keV. For 
comparison, the red region represents a step function at 
the predicted Seitz threshold, where the range represents 
the lcr band including the thermodynamic uncertainty 
and the scale uncertainties in the absolute energy scale 
of the experiment. 



FIG. 3. (Color online) The data points represent the mea- 
sured iodine nucleation efficiency as a function of iodine equiv- 
alent recoil energy, where the contributions from carbon, fluo- 
rine and inelastic scatters have been subtracted. The gradual 
turn on is predominantly due to the angular resolution of 
the experiment, as illustrated by both the red region, rep- 
resenting the step function model with the threshold var- 
ied within the uncertainty on the Seitz theory prediction, 
Et = (13.6 ± 0.6) keV ± 6%, and the dashed blue curve, rep- 
resenting the best fit step function with Et = 16.8 keV. The 
black curve shows the best fit exponential model with the 
threshold constrained by the theory as described in the text. 
The inset shows 2a contours for a fit to the exponential model 
with the threshold allowed to float (pink) or constrained by 
the theory (solid cyan). The colored dots represent the cor- 
responding curves in the main plot. 



Given the energy resolution smearing induced by MCS 
in this experiment, the preference for a value of Et 
higher than the prediction cannot be easily distinguished 
from an exponential model like Eq. (rol) for iodine nucle- 
ation efficiency with a lower threshold energy and a finite 
value of a/ . Previous studies have shown that the Seitz 
model accurately predicts the threshold at which bub- 
ble nucleation begins for heavy radon daughter nuclei in 
CF 3 I [H [10] . We therefore perform a second fit apply- 
ing the exponential model to iodine recoils, taking the 
Seitz threshold calculation as an external input to the 
analysis to explore the allowed range of a/. The best fit 
is shown as the black curve in Fig. [3] The inset shows 
2ct contours for fits to the exponential model with the 
threshold constrained by our prediction (shaded region) 
and free (unshaded region), along with the best fit points. 

To assess the systematic errors associated with carbon 
and fluorine recoils, we refit the data with two alternative 
models for carbon and fluorine efficiency: the flat model 
from jlOj with energy- independent nucleation efficiency, 
VC,F — 49%, above threshold, and a step function with 
Vc,f = 100%. The latter case represents the worst pos- 
sible scenario for the response of the bubble chamber to 



iodine recoils, as the response to carbon and fluorine is 
maximized. We use the exponential model for iodine re- 
coils, allowing a/ to float and treating Et as a nuisance 
parameter constrained by the prediction. The results of 
these fits are summarized in Table [TJ Extended fits over 
the energy interval 9 < E] e < 100 keV have a negligible 
effect on the iodine fit parameters but disfavor the flat 
C/F efficiency model with r\c,F = 49%. 
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TABLE I. Summary of fits to oii, including 90% lower limits 
on Q/. The three different C/F efficiency models described 
in the text are tested, and in all cases the predicted Seitz 
threshold is treated as a nuisance parameter. By maximiz- 
ing the subtracted C/F contribution, the step function with 
tjc,f — 100% represents the worst case for iodine efficiency. 

In conclusion, we have directly measured the efficiency 
for iodine bubble nucleation in a CF3I bubble chamber 
operated with a nominal threshold of (13.6 ± 0.6) keV. 
For some models of carbon and fluorine efficiency, the re- 
sponse to iodine recoils is consistent with a step function 
at the Seitz threshold, but in all cases there is a prefer- 
ence for either a softer turn on or a slightly higher thresh- 
old. Even in the worst case scenario for iodine, however, 
the response of the chamber to iodine recoils is much 
closer to the nominal Seitz model than it is for carbon 
and fluorine recoils. This was expected from the consider- 
ably larger stopping power of iodine, which facilitates the 
concentration of energy that leads to critical bubble for- 
mation. Systematic uncertainties from both the absolute 
beam momentum calibration and the carbon and fluorine 
response limit the present measurement. This measure- 
ment provides confirmation of the sensitivity of COUPP 
bubble chambers to spin-independent WIMP interactions 



with iodine nuclei, a confirmation that was not attainable 
using standard neutron source calibrations. The tech- 
nique of employing hadron elastic scattering as a tool to 
measure bubble nucleation thresholds is now established, 
enabling the measurement of bubble nucleation energies 
on an event by event basis. We have begun studies of the 
feasibility to repeat this technique with different fluids. 
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